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Abstract 

The model at finite temperature is simulated on the lattice. 
For fixed Nt we compute the transition line for N s — » oo by means 
of Finite Size Scaling techniques. The crossings of a Renormalization 
Group trajectory with the transition lines of increasing Nt give a well 
defined limit for the critical temperature in the continuum. By con- 
sidering different RG trajectories, we compute T c /g as a function of 
the renormalized parameters. 



DFTUZ/96/21 DSF 96/50 |hep-lat/961003S 



1 Introduction 



Finite temperature phase transitions in Relativistic Quantum Field Theories 
have been the subject of an intense study since the time of the pioneering 
work of Kirzhnitz and Linde Weinberg and Dolan and Jackiw [[J, 
which revealed that, in general, spontaneously broken symmetries get re- 
stored above a certain critical temperature T c Q While the nature of this 
phenomenon is by now well understood, an accurate description of the tran- 
sition, including sometimes its order, as well as a precise determination of 
the critical temperature can be difficult because of the appearance, at the 
transition point, of infrared divergencies which cause the breakdown of or- 
dinary perturbation theory ||. Even though new analytic non-perturbative 
methods, allegedly immune of infrared problems, have recently been devel- 
oped and have been used to attack this sort of problems |J, the failure of 
the standard analytic methods has motivated some authors to study these 
phase transitions on computers and these days massive simulations of the 
electroweak phase transition are being performed (see for example [[^). 

In this paper we report on the results of a lattice Monte Carlo simulation 
of the A03 theory at finite temperature. There are definite reasons to be 
interested in this simple system. First of all it is believed that the theory 
admits a non-trivial continuum limit (at T = 0) || |9], [U], IT[ [12|], contrary 



to <p\. Moreover, <p\ is superrenormalizable |13| , so that primitive ultraviolet 
divergences occur only in a finite number of diagrams, and this allows one to 
exactly compute its f3 functions. The model is expected to undergo a second 
order phase transition at a finite critical temperature but despite its simple 
ultraviolet properties, its behavior near the transition is particularly hard 
to explore, because infrared divergencies get worsened in three dimensions 
and as a matter of fact there is no reliable way of computing the critical 
temperature in perturbation theory, not even to lowest order WM. Finite 



Temperature Renormalization Group methods may nevertheless be used to 
get estimates of the critical temperature JTH| . 

It appears natural to see if better results can be obtained from the lattice. 
In this respect, there exist theorems proving that on the lattice spontaneously 



1 There exist models, however, which do not seem to follow the general rule and so 
exhibit symmetry breaking at arbitrarily large temperatures |2|, [I). Whether this is a 
genuine effect or an artifact of (lowest order) perturbation theory is still a subject of 
discussion § || @. 
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broken symmetries get always restored at sufficiently high temperatures (in 



an arbitrary number of dimensions) |16|, 17 1. Unfortunatly, these theorems 
only provide upper bounds on the critical temperature (for finite lattice spac- 
ings), which cannot be used to determine the critical temperature in the 
continuum limit, not being strong enough to even guarantee that the critical 
temperature remains finite when the lattice spacing is taken to zero. 

Of the existing lattice studies of the ^-theory at finite temperature, some 
combine approximate analytical treatments with numerical work to get es- 



timates of the critical temperature in the continuum limit Q18], |19| , while 



others use a Monte Carlo simulation pO] . This paper represents a significant 



improvement and extension of these results. We have carried out a Monte 
Carlo simulation on large lattices with massive statistics, making a detailed 
study of finite size effects and of the continuum limit. 

Within the path integral formalism, finite temperature is introduced by 
compactifying the euclidean time direction, temperature being the inverse of 



the temporal extension When doing simulations on the lattice, this goal 
is achieved by using asymmetric lattices N t x N s 2 shorter in the temporal 
direction. However, in order to avoid finite size effects and approach the 
continuum limit, these steps should be followed: 

1. For a finite Nt take the limit N s — > oo. 

2. Take the lattice spacing a to and N t to oo. 

Our aim is to measure the critical temperature for which symmetry 
restoration arises. For fixed N t we have computed, for several values of 
N s , the transition line that separates the broken and symmetric phases. By 
means of Finite Size Scaling (FSS) techniques we have obtained the line for 
N s — > oo. We have done this for several values of Nt. When increasing Nt, 
the transition lines approach that of zero temperature. These critical lines 
have been used to measure T c in the following way. A possible continuum 
theory is defined by a "curve of constant physics" (CCP): it is parametrised 
by the lattice spacing a and, as a — > 0, it approaches the critical line of 
the symmetric lattice. While doing so, it intersects the critical lines of the 
lattices with finite N t . Each intersection point gives a critical temperature, 
T Nt = and T c is just the limit (if it exists) of all these temperatures. 

Due to the super-renormalizability of the model, the asymptotic form 
of the CCP's that approach the gaussian fixed point can be computed ex- 
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actly in perturbation theory. We have considered several CCP's and have 
accepted only those for which the correlation length, measured numerically, 
has shown a satisfactory scaling. For these trajectories we have found a well 
defined limit for the critical temperatures and we have been able to verify the 
relation predicted in |15] between T c and the renormalized parameters of the 
continuum model (defined using dimensional regularization in the Minimal 
Subtraction (MS) scheme) 

T-logf^=- W ^ (TC) , (1) 
V 9 J 9 

where itims{T c ) is the value of the running mass at the scale T c , g is the 
coupling constant and k is a factor not calculable perturbatively. We have 
obtained an estimate of k = 2.1(2). 

The plan of the paper is as follows. In Sec. 2 we present the model. Sec. 3 
defines the curves of constant physics. In Sec. 4 we discuss the effects of 
a finite temperature and explain how to measure its critical value. Sec. 5 
describes the details of the Monte Carlo simulation and the observables that 
have been measured. The results are presented in Sec. 6. Finally, Sec. 7 is 
devoted to the conclusions. 



2 The \cj)\ model 



We consider the theory for a real scalar field in 3 euclidean dimensions, 
described by the bare (euclidean) action: 

, ,o 1 



S 



d A x 



-(^^ + -m^i + f^ 



(2) 



When regularized on an infinite three-dimensional cubic lattice of points 
Q with lattice spacing a, the above action is replaced by the discretized 
expression: 



xen I z ^ L 4 - 

where Aj^ is the lattice derivative operator in the direction /x 



(3) 



<f> L (x + jia) - $ L (x) 
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It is convenient to measure all quantities in the action (^) in units of the 
lattice spacing. We thus define: 



4>l(x) = a 1 ^ 2 $L(x), u = ags, r = a 2 m B . (4) 
In terms of the rescaled quantities the action now reads: 

^ = E{^E(^WW + ^l(x) + |0l(a;)| . (5) 

In order to conform with the customary lattice notation, we perform a 
further redefinition of the field and parameters in ([5]) : 

, . , / _, 24A 1 - 2A - 3/t 
<Pl{x) = V K( Pr, u=—, r = 2 (6) 



and thus write the action in our final form: 

reZ 3 I M J 

This model has been intensely studied at zero temperature, using renor- 
malization group and high temperature expansion techniques [[Kj], and also 



numerically 0, |9], |TT], I2]| , and its features are well known. As it is apparent 
from the action, it has a Z2 symmetry which can be spontaneously broken, if 
the double-well potential is deep enough. In the parameter space (k, A) there 
exist two fixed points: one is the gaussian fixed point (k = 1/3, A = 0) while 
the other is in the same universality class as that of the 3d Ising model. The 
gaussian fixed point has a null attraction domain and the whole transition 
line (except just the origin) possesses the Ising critical exponents. 



3 The Curves of Constant Physics 

In this section we define the CCP's which define the continuum limit of our 
lattice model. They are Renormalization Group trajectories of the theory at 
zero temperature. 

Suppose we start from a symmetric A^ 3 lattice with spacing ao and bare 
parameters m B , gs- In order to take the continuum limit, we have to vary a 
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in the interval [ao, 0) and simultaneously adjust the values of m 2 B , gs, in such 
a way that physical observables approach a constant value when a — > 0. If 
there were no need for renormalization, this would be a trivial task: it would 
suffice to hold m 2 B and gs constant in eq.(f|) and this would tell us how to 
adjust the dimensionless lattice parameters u and r with a in the approach to 
the continuum. Unfortunatly, things are not so easy and the above procedure 
would fail, for example, to produce in the continuum a finite physical mass 
for the lightest boson. Because of the interaction, it is necessary to give the 
lattice parameters a non trivial dependence on a in order for the observables, 
like the masses, to approach a finite limit for a — > and we call "curves of 
constant physics" (CCP's) any dependence of the lattice parameters (w,r) 
(or equivalently (k, A)) on a which achieves this goal. 

In principle, a possible way of defining the CCP's would be by requiring 
that two observables (0\, C 2 ), constructed out of the Green's functions of the 
theory (for example the physical mass of the boson and the coupling constant 
at some fixed momentum scale), keep a constant value in physical units when 
a is varied. The CCP will depend of course on the choice of (Ci,^), but 
for a sufficiently small they will all coincide and define the same continuum 
physics, possibly parametrized in different manners. 

As from a numerical point of view this strategy is very hard to implement, 
we have utilized a different procedure to construct our CCP, based on the fact 
that our model is perturbatively super-renormalizable. Due to this property, 
the perturbative expansion of the theory contains only a finite number of 



ultraviolet primitively divergent diagrams [13]. Concretely, depending on 
the regularization scheme, such divergencies can occur only in the one-loop 
"tadpole" diagram and in the two-loop "sunset" diagram, and both can be 
reabsorbed via a redefinition of the mass term. As a consequence, no infinite 
renormalization of the field and of the coupling constant are needed and we 
can thus identify the bare field and the bare coupling g B in (|J) with the 
renormalized ones, and g^, respectively. We can then rewrite the action 
m as: 



S = d 3 x 



(8) 



_2 y ' 2 a n 4! 
(From now on we shall omit writing the subscript R of all renormalised 
quantities, for simplicity.) In this expression, the only divergent quantity is 
thus m 2 B . We split it as: 

m 2 B = m 2 + 5m 2 (9) 
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where m 2 represents the tree level mass term and Sm 2 is a mass-counterterm 
which diverges in the limit a — > 0. The expression of 8m 2 depends of course 
on the renormalization scheme. Its computation becomes especially easy in 
the lattice analogue of the minimal subtraction scheme, because all is needed 
then are the divergent parts of the lattice tadpole and sunset diagrams. All 
the formulae needed for this computation can be found, for example, in J7| 
and here we just quote the answer (for the limit N —>■ oo): 

m\ = m 2 (y) - — g - ^r^g 2 \n(ay) , (10) 

where y is an arbitrary mass scale and S = 3.1759114. 

With the above expression for the bare mass, the perturbative expansion 
of the Green's functions of the theory is finite to all orders and thus all the 
observables that can be constructed out of them have a well defined limit for 
a — > 0. This suggests that we define as CCP those that simply have a fixed 
value of m 2 (at some fixed scale y) and of g. A similar approach was used 
and shown to work properly for 2d A0 4 f52| . 

Using eq. (|II]) and the definition of the lattice parameters eq.([|), it is 
easy to see that the CCP, starting from some initial point do, To (at the scale 
yo), uq, takes then the form 

u(s) = SUq , 

r(s) = s 2 (r„ + |U (l " j) " 9^2"o lns ) . (") 

where s = a/a Q . Notice that, according to these formulae, the CCP's ap- 
proach the gaussian fixed point, as it must be having been computed within 
perturbation theory. 

The parameters m 2 and g just represent renormalised quantities and do 
not constitute a pair of observables that can be measured directly on the 
lattice. Consider instead one such observable, for example the physical mass 
of the lightest particle, m p h ys , which is simply related to the correlation 
length at large separations £ p % s , since m v h ys = £,ph ys - Now, when computed in 
perturbation theory starting from (g), with m 2 B given by fllUp, m 2 hys becomes 
a function of m 2 , g and a, admitting a finite limit for a — > 0. On dimensional 
grounds we can write this function as 

m lhys = m 2 F(g/m, am). 
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When the continuum limit is sufficiently near, the dependence on a of 
the r.h.s. becomes negligible, m p h ys and thus the correlation length ^ p h ys 
become independent on a and this implies that the lattice correlation length 
£ — iphys/O' becomes proportional to 1/a. 

In terms of the parameter s introduced earlier, the region of continuum 
physics, along the CCP, is that for which the lattice correlation length scales 
according to 

C{s) oc ~ . (12) 

Moreover, one should always have £ > 1 in order for the lattice details to 
be irrelevant and £ < iV in order for finite size effects to be negligible. 



4 Finite temperature 

We briefly discuss the effects of a finite temperature on our system and how 
to measure T c on the lattice. 

Let us suppose that the system is in the ordered phase at T = and 
imagine heating it up. The intuitive picture, valid in general, is that disorder 
increases, as the temperature is raised, until one reaches a critical value 
T c above which no order is possible any more and the full symmetry of 
the lagrangian is restored. This simple picture is confirmed by analytical 
treatments in four space time dimensions [@, || . 

In three dimensions things are less simple. Consider our 0§ model: while 
the existence of the symmetry restoring phase transition at a finite temper- 
ature seems sure, a reliable perturbative computation of the critical tem- 
perature, using standard techniques, is problematic due to the presence of 
severe infrared divergencies. For example, when taking a high temperature 
expansion of the one-loop effective potential one finds that the leading field- 
dependent term becomes complex for small fields, and this leads to a physi- 
cally unacceptable complex value for the critical temperature. The situation 
does not improve upon resumming the so called daisy and super-daisy dia- 
grams [H] : the resulting gap equation determining the thermal mass becomes 



meaningless at the critical temperature [14| . 

To our knowledge, the best analytical estimate of T c for the 0| model is 
the one recently obtained in [|15j using renormalization group methods. It is 
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given in terms of an implicit equation for T c : 



\ 9 

where m 2 MS (fj,) is the running mass 



glogf^U- 8 ™^, (13) 



m 2 MS {n) = m 2 MS (iii) + log ( j • (14) 



According to |L5|, the uncertainty on the value of T c due to the breakdown 
of perturbation theory in the vicinity of the phase transition only affects the 
value of the number k in eq.(|T3"D. It cannot be computed analytically but is 
expected to be of order one and, as will be shown later, we have been able 
to measure it quite accurately. 

The difficulties on the analytical side motivated us to study the phase 
transition of <p\ on the lattice. In order to simulate a finite physical tem- 
perature we work on an N t x iV s 2 lattice. For every N t considered, we have 
performed a series of simulations for several values of N s , and using FSS 
techniques we have obtained the transition line in the thermodynamic limit 
N t x oo 2 . Repeating this analysis for different values of N t we have obtained 
several transition lines. For decreasing N t the critical lines shift deeper inside 
the broken region of the symmetric lattice, though all of them meet at the 
gaussian fixed point (k = 1/3, A = 0) whenever N t > 1. 

Consider now a CCP, call it 7, lying in the ordered region of the T = 
theory and let P be its starting point. Since the critical lines for finite N t 
approach, in the limit Nt — > 00, that of the T = theory, P will also lie in 
the ordered regions of the lattices with finite Nt, for Nt large enough (how 
large it depends on P, of course). Choose now one such Nt, and imagine 
moving along the CCP. In this process a diminishes and thus the physical 
temperature associated with the N t chosen grows up since T = 1/N t a. At 
some point our CCP will cross the transition line for this value of N t and 
symmetry will be restored. If ajv t is the value of the lattice spacing along the 
CCP at the crossing point, we have an estimate of the critical temperature: 

T Nt = j^— (15) 
N t a Nt 

In order to get the continuum critical temperature we have to let N t — > 00. 
When Nt is increased, the corresponding critical line gets closer to that of the 
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T = theory and 7 will intersect it for a smaller value of a. The temperatures 
T/v t are thus expected to approach a limit, which is the desired continuum 
critical temperature T c for the continuum theory associated with 7. 

In the absence of an external input which fixes the energy scale, we can 
assign a value a = a to an arbitrary point on the CCP and then the value of 
a at any other point on the trajectory is given by equations (pp . In this way 
we can compare the temperatures T Nt corresponding to different values of N t 
for a given CCP, but it is not possible to put together the values from distinct 
CCP's due to the arbitrariness in the scale. To avoid this problem, one can 
measure adimensional quantities, obtained for instance from the quotient of 
two observables with the same dimensions. Of this sort are the quantities 
T c /g and m 2 MS {T c ) / g 2 entering equation (|Hf). 

The analysis of jn| was carried in the continuum, using dimensional reg- 
ularization and minimal subtraction, while we are using the lattice regular- 
ization. So, in order to compare the results of our simulations with eq. ([13|) , 
we have to first relate the continuum renormalized parameters appearing in 
eq.(|T3|) to ours. In view of what we said earlier about the ultraviolet proper- 
ties of the model, the only non trivial relation that needs to be computed is 
the one among the mass parameters of the two schemes. Using the formulae 
of for the relevant two-loops diagrams contributing to the self-energy, it 
is not hard to check that: 

m 2 M M = m 2 (/x) - ^(log6 + C) + o(g 3 ) , (16) 

where ( = 0.09 and m 2 (/i) is the lattice renormalized mass defined in eq. ([T0|) . 
Upon using egs-fliTf), (PH) and ([To]), we can now relate the adimensional 
quantities T c /g and m 2 MS (T c ) / g 2 to the simulation parameters according to 
the following equations: 

= Jim -Ar + ^u- ^- 2 u\\o g N t + log6 + C)) , (17) 

rpc I 

— = lim (18) 

g Nt^oo N t U 

where u and r are evaluated at the points in which the CCP intercepts the 
critical line of the iV r lattice. In our simulation we have measured the above 
quantities for various CCP's, and this allowed us to check the validity of 
eq. ([l3|) and measure, at the same time, the value of the unknown parameter 
k. 
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5 Monte Carlo simulation and observables 



In the course of the simulations, most of our efforts went into the computation 
of the critical lines for different values of N t . On the other hand, we have 
made measurements also on the symmetric lattice in order to identify the 
scaling region of our perturbative CCP's. 

We have used a Metropolis algorithm combined with a Wolff single clus- 
ter method |23| (the latter updates the sign of the field), taking a ratio of 
20 clusters every 3 Metropolis iterations or 10 clusters every 2 Metropolis. 
These values are inside the region of ratios we have found nearly optimal for 
decorrelation versus CPU time close to the critical lines. When, in the intent 
of checking the CCP's, we have measured the correlation lengths on the sym- 
metric lattices, the simulations have been performed in the ordered region of 
the phase diagram and there the clustering is not so effective. However these 
simulation points did not lie very deep in the broken phase and then we have 
used this method as well. 

The most intensive simulations have been those relative to the asymmetric 
lattices. We have used four different values of N t (2,3,4,6) and for each of 
them we have computed the critical value of k for four values of A (0.0005, 
0.001, 0.002, 0.003). For each pair (N t ,X), we have considered lattices of 
N s = 12, 16, 24, 32, performing on each of them a single simulation for one 
value of k. This value has been chosen near the critical one, by means of 
an hysteresis or by extrapolating from the values found for the previous 
lattices. Afterwards, we have used the spectral density method (SDM) |24f 
to extrapolate the measurements to other values of k in the vicinity of the 
simulation point. Finally, the value of K-critical in the limit N s — > oo has 
been computed by examining the crossings of the Binder cumulants relative 
to lattices with different values of N s , as explained below. 

At each of the points and lattices that have been simulated, the number 
of iterations made has been approximately 1, 500, 000 x (3 Metropolis + 
20 Clusters) or 2, 500, 000 x (2 Metropolis + 10 Clusters), for a total CPU 
time equivalent to nearly one year of workstation. Most of the simulations 
have been carried out with iterations consisting of 3 Metropolis+20 Clusters 
and performing measurements only at intervals of 10 iterations. The errors 
in the estimation of the observables have been calculated with the jackknife 
method. 

We have measured the integrated autocorrelation times r. The values ob- 
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Figure 1: Logarithms of the autocorrelation times r as a function of 
log(iVtiV fl 2 ) 1/3 , for N t = 3 (full circles), N t = 4 (empty triangles) and iV t = 6 
(full squares) and N s = 12, 16, 24, 32. The straight lines are the fits for the 
dynamical exponents z according to eq.flT§p. The points not covered by the 
fits correspond to N s = 12 and have not been used for the computation of 
the fits. 

tained from the K-energy, A-energy and magnetisation are totally compatible. 
They oscillate, in units of 10 x (3 Metropolis + 20 Clusters), approximately 
from r ~ 3 in the best case (largest A and smallest lattice) to r ~ 38 in the 
worst (smallest A and largest lattice). In particular we have estimated the 
dynamical exponent z defined by 

t ~ (N t x N s 2 f 3 (19) 

with fixed Nt, obtaining z = 0.8(1). This estimate is not totally rigorous be- 
cause the autocorrelation times have been measured at the simulation points 
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which do not coincide exactly with the 'critical' ones (but they do it approx- 
imately). However it should be basically correct because, as can be seen in 
Figure [TJ, a power law with a rather well defined exponent is observed for the 
largest lattices. Figure 1 shows how r scales with N s , for fixed N t = 3,4,6 
at A = 0.001 and A = 0.003. The fits have been computed using only the 
data relative to N s = 16, 24, 32, but for completeness in the figure we have 
displayed also the points corresponding to N s = 12. They give an exponent 
z oscillating between 0.7 and 0.9 approximately. Metropolis alone would give 
a value z = 2 while the Wolff algorithm, in spin systems where can it be used 
by itself, has z = because at the critical point the clusters can have any size. 
In our case the clustering is only able to change the sign of the fields, but its 
non-locality contributes to diminish significantly the dynamical Metropolis 
exponent. 

For the sake of completeness, we have also determined the approximate 
position of the critical line at zero temperature. In this case no extrapolation 
to the thermodynamical limit has been made, and we have simulated only a 
16 3 lattice. 

In order to estimate the critical points, we have used the Binder cumulant 
defined as: 

u "-w = h$Fr (20) 

where 

M=\{l/V)^r\ (21) 
r 

is the magnetisation or expectation value of the field, V = N t x N% being 
the volume. 

At several points on the CCP's, we have measured the second moment 
correlation length on the symmetric lattice 16 3 , in order to check the scaling. 
It is defined as |H l26l 




e= hr r^ , (22) 



where G(x) is the correlation function. Usually this quantity is measured in 
the disordered phase, while we need evaluate it in the broken phase. There 
(M) and then we have to use the connected correlation function 
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G(x,y) = G(x -y) = (<f> x <j> y ) - (<j> x ) {<j> y ) . (23) 
Using the lattice symmetries, £ can be expressed as 

where N = N s = N t is the lattice size, x is the susceptibility 

X = E^)=n(M 2 )-(M) 2 ) (25) 

X 

and F is the analogous quantity at the smallest nonzero momentum (having 
two null components and one equal to ±2n/N) 



F = Y, G(x)e ip - X = 7^(1 E <Pre 2mri/L \ 2 + | ]T r e 2 ™ r2 / L | 2 + | £ <p r e 2mr * ,L \ 2 ) 

X V V V 

(26) 

r-j being the components of r (indeed, the fact of being in the ordered phase 
is irrelevant as far as F is concerned, because the disconnected part cancels 
out from the sum due to the Fourier factor). 



6 Results 

6.1 Phase diagram 

As explained in the previous section, for each pair (N t , A) we have used 
iV s = 12, 16, 24, 32. As the field is not compact, it is difficult to obtain the 
critical points from observables like the specific heat and the susceptibility. 
An estimate of the critical point for each value of N s could be obtained from 
the maximum of the derivative of the cumulant Un s (k). However, in order to 
minimize finite size effects, we have rather used the intersections k*(N Si , N S2 ) 
among all possible pairs of curves U Ns (k,), because the differences between 
them come only from the corrections to scaling. After simulating one value of 
k for every N s (having fixed Nt and A), we have computed U n s (k) in a narrow 
interval around the simulation point by means of the SDM. The simulation 
point was known to be near the critical line thanks to some hysteresis or from 
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Figure 2: The Binder cumulants Un s as a function of the hopping parameter 
k, for the lattices with N s = 12 (circles), N s = 16 (triangles), N s = 24 
(squares), N s = 32 (diamonds). The full marks are the simulation points, 
while the empty ones are extrapolations via the SDM. 



the estimations for the lattices simulated previously. Figure || shows Un s («) 
for the different values of N s , with N t = 2 and A = 0.003, as an example. 
The full marks correspond to the simulation points, while the empty ones 
are the extrapolated values. 

Starting from the intersections k*(N Si , N S2 ), the critical value k c in the 
thermodynamic limit N s — » oo can be estimated by using the following law 



k*(N s , bN s ) - k c = LA^jV,-"- 1 ^ (27) 

where u is the exponent for the corrections-to-scaling. Since the scaling 
parameter is iV s , while Nt is fixed, according to the hypothesis of dimensional 
reduction and universality, we have used the exponents of the Ising model in 
two dimensions, namely u — 1, a>~4/3 fT5| , obtaining satisfactory fits. To 
give an idea of their quality, in Figure |3| we show two extreme cases, smallest 
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Figure 3: Intersections of Binder cumulants K*(N s ,bN s ) as a function of 
\Z/u" 1 Nj u ~ x ' v (squares). The dotted lines are the fits according to eq. (|27|). 
The full circles are the extrapolated values for k c in the limit N s — > oo. 



A - largest N t and largest A - smallest N t , which are among the worst and 
the best fits respectively. 

On the other hand, as a further justification for the use of the two- 
dimensional exponents, we have studied the scaling of the maximum deriva- 
tive of the Binder cumulant with N s . It was expected to scale like N^ u and 
we have checked this behavior for every pair (N t , A), obtaining a good scaling 
and values of v close to 1 in every case [v ranges from 0.95 to 0.99, and is in 
all cases almost compatible with 1 within the errors). Figure [| shows such a 
fit for N t = 3, A = 0.001. 

Although we have not made an exhaustive study of the order of the tran- 
sitions, all the indications point towards the expected second order. One 
element is the value of the exponent v obtained. Another is that no clear 
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log(N) 

Figure 4: Logarithms of the maximal Binder cumulant derivatives as a func- 
tion of N s for N t = 3, A = 0.001. The line is the fit to compute v. 

discontinuity is observed in the histograms for the lattices that we have sim- 
ulated. In Figure ^ we plot the histograms of the K-energy for N t = 6, 
N s = 16, 24, 32 at the point A = 0.003, k = 0.3364, which is very close to the 
transition line. 

In table p] we quote the values of k c for the different (N t , A) studied. In 
order to have a continuous transition line for every value of N t , we have made 



A 


Kc (N t = 6) 


*c{N t = 4) 


K c (N t = 3) 


K c (N t = 2) 


0.003 


0.336211(30) 


0.336648(17) 


0.337024(30) 


0.338170(29) 


0.002 


0.335249(19) 


0.335567(27) 


0.336007(27) 


0.336711(24) 


0.001 


0.334373(9) 


0.334597(21) 


0.334751(26) 


0.335243(45) 


0.0005 


0.333835(23) 


0.333924(19) 


0.334024(35) 


0.334237(16) 



Table 1: Critical values k c for the different values of N t and A. 
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Figure 5: Hystogram of the configurations distribution as a function of 
the term in the action eq.(^) proportional to the hopping parameter k, for 
N t = 6, A = 0.003 and k = 0.3364, N s = 16 (dots), 24 (dashes), 32 (solid). 

a linear extrapolation between each pair of correlative A values. Adding the 
critical line computed for the symmetric lattice (estimated on a 16 3 lattice 
at A = 0.001,0.003 and then continued up to the gaussian point) we have 
obtained the phase diagram shown in Figure |6], where we have also included 
the CCP's we have used for the determination of T c . 

6.2 Critical temperatures 

Now we have to search for the intersections of some CCP's with the transition 
lines for the different values of N t . The region of the phase diagram that we 
have computed limits the window of accessible CCP's, as we could only use 
CCP's whose intersections with the critical lines are in the range 0.0005 < 
A < 0.003. In addition, in order to ensure that we were exploring a region 
of continuum physics, we have required a good scaling of the correlation 
length along the trajectories. Trajectories below D in Figure |6| intersect only 
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Figure 6: Critical lines (solid curves) for N t = 2,3,4,6. The dashed line 
is the 'critical' line for the symmetric lattice 16 3 . The dotted lines are the 
CCP's A,B,C,D. 



the critical lines with Nt = 2,3 and then have been discarded. Trajectories 
above A have shown a poor scaling and have also been eliminated. This 
could depend on two circumstances: on one side the scaling region might be 
too near the gaussian point for these trajectories, on the other it is quite 
possible that near the critical line the finite size effects in the measurement 
of the correlation length are important. 

These considerations have led us to consider only the trajectories included 
in the region of the phase diagram bounded by the curves A and D in Figure 
|6|. All in all, we have used four trajectories, A, B, C and D, starting from the 
points with A = 0.003 and « = 0.3382, 0.3389, 0.3395, 0.3402 respectively. 
Among them, the only one cutting the critical line for Nt — 6 is A, but it does 
it almost tangentially and the error in the determination of the intersection 
point is consequently very large. This produces an error still larger in the 
estimate of the temperature because in the proximity of the gaussian point 
the lattice spacing a changes very rapidly along the CCP. As a matter of 
fact we have not been able to use Nt — 6 for the measurement of the critical 
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AT 


ft 


A 


? 


? • s 


-1671-^(7^) 


A7rT Nt 

a 


2 


0.33810(8) 


0.00296(3) 


3.42(3) 


3.38(3) 


17.275(2) 


10.11(11) 


3 


0.33584(9) 


0.00187(5) 


5.78(9) 


3.65(6) 


17.268(4) 


10.55(28) 


4 


0.33489(7) 


0.00130(4) 


8.84(12) 


3.91(6) 


17.256(5) 


11.29(35) 


2 


0.33722(6) 


0.00236(3) 


3.75(5) 


2.98(4) 


22.389(2) 


12.64(14) 


3 


0.33536(8) 


0.00148(5) 


7.10(14) 


3.59(7) 


22.381(5) 


13.23(38) 


4 


0.33469(6) 


0.00110(3) 


9.36(16) 


3.53(6) 


22.381(5) 


13.29(39) 


2 


0.33674(6) 


0.00202(3) 


3.94(8) 


2.70(6) 


26.829(2) 


14.68(17) 


3 


0.33508(8) 


0.00127(4) 


7.27(32) 


3.15(15) 


26.821(5) 


15.45(48) 


4 


0.33456(9) 


0.00097(5) 


9.92(28) 


3.31(10) 


26.825(9) 


15.06(86) 


2 


0.33634(7) 


0.00175(3) 


4.11(8) 


2.45(5) 


32.074(3) 


16.95(26) 


3 


0.33486(7) 


0.00109(4) 


8.02(25) 


3.02(10) 


32.065(6) 


17.91(60) 


4 


0.33427(12) 


0.00076(7) 


10.93(19) 


2.86(6) 


32.052(16) 


19.28(1.89) 



Table 2: Correlation lengths £ and estimates T^ t of the critical temperatures 
associated with the crossings of the CCP's A, B, C, D (from top to bottom) 
with the critical lines of the lattices with N t = 2, 3, 4. 

temperature and we had to limit ourselves to the crossings of the CCP's with 
the critical lines for Nt = 2, 3, 4. 

In table ^| we show the values of k and A at the points where the tra- 
jectories A, B, C and D (from up to down) intersect the critical lines for 
N t = 2,3,4. Their errors come from the deviations in the estimation of the 
critical lines, and have been calculated as explained below for the critical 
temperatures. 

In correspondence with these intersection points we have measured the 
correlation length on a 16 3 lattice in order to verify its scaling. We have 
chosen these and not any other points on the CCP's because the good or 
poor scaling of £ just there indicates which values of N t we will be able to 
use for estimating the critical temperature in the continuum. The values of 
£ obtained are quoted in table [| The errors quoted are the ones arising from 
the simulation and have been calculated with the jackknife method. 

When going towards the gaussian point, £ should grow inversely to the 
parameter s, i.e. inversely to a (see. eq. ([l2]) ) and thus the product £ • s 
should be constant along the trajectory (the absolute scale of s is of course 
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arbitrary and, for each CCP, we have taken s = 1 at its starting point). It 
is apparent from table § that for the trajectories B, C and D the values of 
£ • s for N t = 3, 4 are fully compatible within the errors, while for N t = 2 
there are significant deviations from scaling. For the trajectory A the result 
is less satisfactory, something that will be reflected also in a small difference 
among the estimates for the critical temperature resulting from N t = 3 and 
4 for this trajectory. 

As we said above, we have constructed the critical lines by means of a 
linear interpolation of the critical points obtained numerically and have then 
used them to estimate the intersections with the CCP's. From Figure ^| it is 
apparent that this is a rather good approximation. We have estimated the 
errors in the measurement of T^Jg and m 2 MS /g 2 in the following way. In 
addition to the line formed by joining the values of k c (X) for the four values 
of A simulated, we have considered two more lines, on the two sides of the 
previous one. One has been obtained by joining the points corresponding to 
the upper values k c (A) + cr(/t c (A)) and the other by joining the points corre- 
sponding to the lower values k c (A) — cr(ft c (A)), o"(k c (A)) being the error in the 
determination of k c (X). Each of the three lines gives a different estimate for 
the observables. We have taken as our estimate the one given by the central 
line and the error has been taken to be one half of the difference between the 
values given by the two external lines. This procedure would overestimate 
the error if the critical lines between the critical points computed were really 
straight; this overestimation should approximately compensate for the error 
introduced by the linearization. 

The values of T^J g and rn 2 MS /g 2 obtained in this way are shown in table 
|2| (multiplied by the constant factors Atc and — lQir 2 in order to compare 
directly with the results in [fL5|| ). It can be seen how, for N t = 3,4, T^Jg 
has reached a well defined limit (only for the trajectory A the values of T^/g 
and T^/g are not perfectly compatible within the errors). 

In order to check eq. (|T3"|) , we have performed a series of best fits. To 
this purpose we have used all the data relative to N t = 3, 4, except the 
point corresponding to the intercept of trajectory A with the critical line 
for N t = 3, for the reason explained above. Some care has been necessary, 
because not all the data are statistically independent. Consider for instance 
the intercepts of trajectories B, C, D with the critical line for N t = 3: they 
belong to the same straight segment of the interpolated critical line, which 
implies that among them only two are statistically independent of each other, 
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while the third can be determined thereof. For this reason, we could not use 
all the data at once in a single fit; rather we have repeated the fits for all 
possible maximal subsets of statistically independent points extracted from 
the data. It turned out that all subsets consisted of five points out of a total 
of seven. The plots of x 2 (^) f° r a U the combinations of data are displayed in 
Figure |7| (notice that x 2 {k) is not normalised per degree of freedom). 

As can be seen from the figure, the errors on the estimates of k (equal 
to semi-widths of the intervals in k, around the minimum, for which the 
variation of x 2 is less than 1) are comparable for all the combinations of data 
and are around 0.15. It is apparent that in all cases the minimum value of 
X 2 {k) per degree of freedom (four) is very good, which indicates that the fits 
can be trusted. The values of k obtained from the fits are all practically 
consistent with each other, within the errors. As our final estimate, we have 
taken their average, while for the error we have taken a value larger than the 



Zi i i i | i i i i | i i i i | i i i i | i i i i | i i i 




1.8 1.9 2.0 2.1 2.2 2.3 2.4 



k 



Figure 7: Plots of x as a function of the parameter k in eq. ([T3|) for all 
maximal subsets of statistically independent data extracted from table |2] 
(the point relative to the CCP A and N t = 3 has been excluded, so as all 
those of N t = 2). The number of degrees of freedom in all cases is 5 — 1 = 4. 
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Figure 8: Estimates T Nt of the critical temperature as a function of the 
renormalised parameters m 2 MS {T]s[ t ) and g, for the CCP's A,B,C, D (from 
left to right), N t — 3 (full circles) and N t = 4 (empty triangles). The solid 
line is the fit according to eq. flT3|) , with k = 2.1. The dotted line corresponds 
to k = 1. 

individual errors, in such a way that all the fits are compatible with it. In 
conclusion, we have got: k = 2.1(2). The final fit is shown in Figure [8],where 
the curve having k — 1 is also displayed for comparison. 

7 Conclusions 

We have investigated the phenomenon of symmetry restoration at finite tem- 
perature in the scalar A03 model, by means of a Monte Carlo simulation. Us- 
ing FSS techniques we have carefully determined the critical lines separating 
the ordered and the unordered phases on asymmetric lattices, for a number 
of different extensions in the temporal direction. The intersections of these 
lines with various RG trajectories lying in the ordered phase of the T = 
theory provide estimates of the critical temperature. The analysis of the 
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results for lattices with increasing extensions in the time-direction clearly 
shows that these temperatures approach a well-defined finite value in the 
continuum limit, confirming the expectation of a symmetry restoring phase 
transition at a finite physical temperature for the continuum theory. We have 
measured T c for various RG trajectories. In this way we have checked the 
validity of the following equation relating T c to the renormalized parameters 
of the continuum theory: 

/47rfcT°\_ 87m4 g (T c ) 
V 9 J 9 
which was derived in |T5] using renormalization group methods, and have 



been able to measure the so far unknown parameter k, obtaining the value 
k = 2.1(2). 
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